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Abstract 

Most of the theoretical physics known today is described by using a small 
number of differential equations. If we study only linear systems, different 
forms of the hypergeometric or the confluent hypergeometric equations of- 
ten suffice to describe this problem. These equations have power series 
solutions with simple relations between consecutive coefficients and can 
be generally represented in terms of simple integral transforms. If the 
problem is nonlinear, one often uses one form of the Painleve equation. 
There are important examples, however, where one has to use more com- 
plicated equations. An example often encountered in quantum mechanics 
is the hydrogen atom in an external electric field, the Stark effect. One 
often bypasses this difficulty by studying this problem using perturba- 
tion methods. If one studies certain problems in astronomy or general 
relativity, encounter with Heun equation is inevitable. This is a general 
equation whose special forms take names as Mathieu, Lame and Coulomb 
spheroidal equations. Here the coefficients in a power series expansions 
do not have two way recursion relations. We have a relation at least 
between three or four different coefficients. A simple integral transform 
solution also is not obtainable. Here I will try to introduce this equa- 
tion whose popularity increased recently, mostly among theoretical physi- 
cists, and give some examples where the result can be expressed in terms 
of solutions of this equation. 



1 Introduction 



Most of the theoretical physics known today is described by using a small num- 
ber of differential equations. If we study only linear systems, different forms of 
the hypergeometric or the confluent hypergeometric equations often suffice to 
describe this problem. These equations have power series solutions with simple 
relations between consecutive coefficients and can be generally represented in 
terms of simple integral transforms. If the problem is nonlinear, one often uses 
one form of the Painleve equation. 

Let us review some well known facts about linear second order differential equa- 
tions. Differential equations are classified according to their singularity structure 
[U [2] ■ If a differential equation has no singularities over the full complex plane, 
it can only be a constant. Singularities are classified as regular or irregular 
singular points. If the coefficient of the first derivative has at most single poles, 
and the coefficient of the term without a derivative has at most double poles 
when the coefficient of the second derivative is unity, this second order differen- 
tial equation has regular singularities. Then we have one regular solution while 
expanding around these singular points. The second solution around a regular 
singular point has a branch cut. If the poles of these coefficients are higher, we 
have irregular singularities and the general solution has essential singularities 
around these points [3]. 

As stated in Morse and Feshbach |4|, an example of a second order differential 
equation with one regular singular point is 



d 



2 



0. (1) 



This equation has one solution which is constant. The second solution blows up 
at infinity. The differential equation 

— + fc^. = 0. (2) 

has one irregular singularity at infinity which gives an essential singularity at 
this point. The equation 

d'^w , , dw , , 

+ + 

has two regular singular points, at zero and at infinity. In physics an often used 
equation is the hypergeometric equation 

z{\ - z)-— -1- c - 1 + a + 5 z — - ahw = 0. (4) 
dz^ dz 

This equation has three regular singular points, at zero, one and infinity. Jacobi, 
Legendre, Gegenbauer, Tchebycheff equations are special forms of this equation. 
When the singular points at z=l and z equals infinity are " coalesced" [S] at 
infinity, we get the confluent hypergeometric equation 

(i^w dw 
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with an essential singularity at infinity and a regular singularly at zero. Bessel, 
Laguerre, Hermite equations can be reduced to this form. 

An important property of all these equations is that they allow infinite series 
solutions about one of their regular singular points and a recursion relation can 
be found between two consecutive coefficients of the series. This fact allows 
us having an idea about the general properties of the solution, as asymptotic 
behaviour at distant points, the radius of convergence of the series, etc. 
A new equation was introduced in 1889 by Karl M. W. L. Heun [5]. This is 
an equation with four regular singular points. This equation is discussed in the 
book edited by Ronveaux [7] . All the general information we give below is taken 
from this book. 

As discussed there [8], any equation with four regular singular points can be 
transformed to the equation given below: 

(Pw c d e dw abz — q 

+ + ^ + w]d7-.(z-l)(z-/)"' = °- 

There is a relation between the constants given asa + b+ l = c + d + e . Then 
this equation has regular singularities at zero, one, / and infinity. If we try to 
obtain a solution in terms of a power series, one can not get a recursion relation 
between two consecutive coefficients. Such a relation exists between at least 
three coefficients. A simple solution as an integral transform also can not be 
found [1]. 

One can obtain different confiuent forms of this equation. When we "coalesce" 
two regular singular points at zero and infinity, we get the confluent Heun 
equation [S]: 

-{{z^ - 1)-) + [-Az^ - 1) + 2pPz - A _ 1) = 0. (7) 

Special forms of this equation are obtained in problems with two Coulombic 
centers [10], 

i - 1)^) + - 1) + 2pPz - A - j^^]^ = 0, (8) 

whose special form is the spheroidal equation |10] . 

A((.^-1,^, + |-pV-1,-A-^|. = 0. (9) 
Another form is the algebraic form of the Mathieu equation [10] : 

i^„,^-i,*) + lV(.'-i)-A-jpL_|„,.o^ (10, 

If we coalesce two regular singular points pairwise, we obtain the double con- 
fluent form [TT] : 

D^w + aiz + ^ + Bi + ^ z + Bo + ^-^ + B-i - ^)-)]w = 0. 
z 2 2 2 z 

(11) 
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Here D — z-^. If we equate the coefRcients ai,i = —1, 1 and set Si, i = — 1, 1 to 
zero, we can reduce the new equation to the Mathieu equation |12j . an equation 
with two irregular singularities at zero and at infinity. Another form is the 
biconfluent form, where three regular singularities are coalesced. The result is 
an equation with a regular singularity at zero and an irregular singularity at 
infinity of higher order }13l : 

^cPw dw , . , . T . c. /I. , . 

^ ^'dl'^ ^ ^ ^ ^ ~ ^ ^' 

The anharmonic equation in three dimensions can be reduced to this equation 
[HI: 

^ + {E-^-^ir^-Xr^^ rjr^)w = 0. (13) 

In the triconfluent case, all regular singular points are "coalesced" at infinity 
which gives the equation below |15j : 

d^w 9 

+ (Ao + Aiz + Asz^ - -z^)w = 0. (14) 

These different forms are used in different physics problems. 

In SCI I found 94 papers when I searched for Heun functions. More than two 
thirds of these papers were published in the last ten years. The rest of the 
papers were published between 1990 and 1999, except a single paper in 1986 on 
a mathematical problem jl6j . This shows that although the Heun equation was 
found in 1889, it was largely neglected in the physics literature until recently. 
Earlier papers on this topic are mostly articles in mathematics journals. 
The list of books on this topic also is not very long. There is a book edited by 
A.Ronveau, which is a collection of papers presented in the "Centennial Work- 
shop on Heun's Equations: Theory and Application. Sept. 3-8 1989, Schloss 
Ringberg" . It was published by the Oxford University Press in 1995 by the title 
Heun's Differential Equations [7|. There is an extensive bibliography at the end 
of this book. I also found frequent references to two other books: Mathieusche 
Funktionen und Sphaeroidfunktionen mit Anwendungen auf physikalische und 
technische Probleme by Joseph Meixner and Friedrich Wilhelm Schaefke, pub- 
lished by Springer Verlag in 1954 [IT and a 1963 Dover reprint of a book 
first published in 1946, Theory and Applications of Mathieu Functions by N.W. 
McLachlan [TH] . These are on functions which are special cases of the Heun 
equation. Some well known papers on different mathematical properties of these 
functions can be found in references [TO] - [53]. 

A reason why more physicists are interested in the Heun equation recently may 
be, perhaps, a demonstration of the fact that we do not have simple problems 
especially in " the theoretical theory" section of particle physics anymore and 
people doing research on this field have to tackle more difficult problems, either 
with more difficult metrics or in higher dimensions. Both of these extensions 
may necessitate the use of the Heun functions among the solutions. We can give 
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the Eguchi-Hanson case as an example. The wave equation in the background of 
the Eguchi-Hanson metric |24| in four dimensions has hypergeometric functions 
as solutions [33] whereas the Nutku helicoid [57] metric, another instanton 
metric which is little bit more complicated than the Eguchi-Hanson one, gives 
us Mathieu function solutions [2H], a member of the Heun function set. We also 
find that the same equation in the background of the Eguchi-Hanson metric, 
trivially extended to five dimensions, gives Heun type solutions [29] . 

Note that the problem need not be very complicated to end up with these 
equations. We encounter Mathieu functions if we consider two dimensional 
problems with elliptical boundaries |30|. Let us use x = ^ a cosh ^ cos 6*, y = 
iasinh/isin^, where a is the distance from the origin to the focal point. Then 
the Helmholtz equation can be written as 

+ deei^ + ^a'^k'^[cosh^ - cos^ 6]^; = (15) 
which separates into two equations 

"^'^ (6 - cos^ e)H = 0, (16) 



(PM 



+ {b-h'^ cosh^ fi)M ^ 0. (17) 



a, fc, h are constants given in j30j . The solutions to these two equations can be 
represented as Mathieu and modified Mathieu functions. 

If we combine different inverse powers of r, starting from first up to the fourth, 
or if we combine the quadratic potentials with inverse even powers of two, four 
and six, we see that the solution of the Schrodinger equation involves Heun 
functions [31j . Solution of the Schrodinger equation to symmetric double Morse 
potentials, like V{x) = ^sinh2x — (s -f ^)Bcoshx , where s = (0, 1/2, 1, ) [35], 
also needs these functions. Similar problems are treated in references [551 [Ml 155] . 

o In atomic physics further problems such as separated double wells, Stark 
effect, hydrogen molecule ion use these functions. Physics problems which end 
up in these equations are given in the book by S.Y. Slavyanov and S. Lay [36] . 
Here we see that even the Stark effect, hydrogen atom in the presence of an 
external electric field, gives rise to this equation. As described in page 166 of 
Slavyanov's book, cited above (original reference is Epstein [37] , also treated by 
S.Yu Slavyanov [55]), when all the relevant constants, namely Planck constant 
over 2tt, electron mass and electron charge are set to unity, the Schrodinger 
equation for the hydrogen atom in a constant electric field of magnitude F in 
the z direction is given by 

(A + 2[£;-(Fz-i)])* = 0. (18) 
Here A is the laplacian operator. 
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Using parabolic coordinates, where the cartesian ones are given in terms of the 
new coordinates by cc = y/^cosfj), y — \/^sincj), z ~ and writing the wave 
function in the product form 

= \/^lV{^)U{r))exp{im<p], (19) 

we get two separated equations: 

fe+^)nC)=0, (20) 









de ^ 


(f- 




d^U 






diff' 







F 1 



-TV+-^mv) = 0. (21) 



Here (3i and /32 are separation constants that must add to one. We note that 
these equations are of the biconfluent Heun form. 

The hydrogen molecule also is treated in reference |36j . (original reference is 
[39] ) . When the hydrogen-molecule ion is studied in the Born-Oppenheimer 
approximation, where the ratio of the electron mass to the proton mass is very 
small, one gets two singly confluent Heun equations if the prolate spheroidal 
coordinates = ^^^c^ , rj = 2liii2. are used. Here c is the distance between the 
two centers. Assuming 



V' = V^vViOU{v)expiim^), (22) 
we achieve separation into two equations: 

|(„_a9 + (A^«'-»«-^ + .)v-.0, (.3, 

i?.(,_,=,«),(,V-^..)a.o. (.4, 

Both equations are of the confluent Heun type. 

If we mention some recent papers on atomic physics with Heun type solutions 
we find three relatively recent papers which treat atoms in magnetic fields: 

o Exact low-lying states of two interacting equally charged particles in a 
magnetic field are studied in by Truong and Bazzali |40| . 

o The energy spectrum of a charged particle on a sphere under a magnetic 
field and Coulomb force is studied by Ralko and Truong [IT]. Both papers 
reduce the problem to a biconfluent Heun equation. 

o B.S. Kandemir presented an analytical analysis of the two-dimensional 
Schrodinger equation for two interacting electrons subjected to a homogeneous 
magnetic field and confined by a two-dimensional external parabolic potential. 
Here also a biconfluent Heun equation is used |42| . 

o Dislocation movement in crystalline materials, quantum diffusion of kinks 
along dislocations are some solid state applications of this equation. The book 
by S.Y. Slavyanov and S. Lay [36j is a general reference on problems solved 
before 2000. 
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o In a relatively recent work P. Dorey, J. Suzuki, R. Tateo [13] show that 
equations in finite lattice systems also reduce to Heun equations. 

In the rest of this paper I will comment only on papers on particle physics and 
general relativity. 

o In a relatively early work, Teukolsky studied the perturbations of the Kerr 
metric and found out that they were described by two coupled singly confluent 
Heun equations [H]. 

o Quasi-normal modes of rotational gravitational singularities were also stud- 
ied by by E.W. Leaver [45_ with Heun type equations. 

In recent applications in general relativity, these equations become indispensible 
when one studies phenomena in higher dimensions, or in different geometries. 
Some references are: 

o D. Batic, H. Schmid, M. Winklmeier studied the Dirac equation in the 
Kerr-Newman metric and static perturbations of the non-extremal Reisner- 
Nordstrom solution and ended up with Heun type solutions [HI. D. Batic and 
H. Schmid also studied the Dirac equation for the Kerr-Newman metric and 
computed its propagator [57]. They found that the equation satisfied is a form 
of a generalized Heun equation described in Reference [IB] . In later work Batic, 
with collaborators, continued studying Heun equations and their generalizations 
MM- 

Prof. P.P. Fiziev studied problems whose solutions are Heun equations exten- 
sively: 

o In a paper published in'gr-qc/0603003' he studied the exact solutions of the 
Regge- Wheeler equation in the Schwarschild black hole interior [SO] . In another 
paper he obtained exact solutions of the same equation and quasi-normal modes 
of compact objects [5T] . 

o He presented a novel derivation of the Teukolsky-Starobinsky identities, 
based on properties of the confluent Heun functions [52] . These functions define 
analytically all exact solutions to the Teukolsky master equation, as well as to 
the Regge- Wheeler and Zerilli ones. 

o In a talk given at 29th Spanish Relativity Meeting (ERE 2006), he depicted 
in more detail the exact solutions of Regge- Wheeler equation, which described 
the axial perturbations of Schwarzschild metric in linear approximation, in the 
Schwarzschild black hole interior and on Kruskal-Szekeres manifold in terms of 
the confluent Heun functions ,^3j. 

o All classes of exact solutions to the Teukolsky master equation were de- 
scribed in terms of confluent Heun functions in Reference [M , . 

o In reference [SS] he reveals important properties of the confluent Heun's 
functions by deriving a set of novel relations for confluent Heun's functions and 
their derivatives of arbitrary order. Specific new subclasses of confluent Heun's 
functions are introduced and studied. A new alternative derivation of confluent 
Heun's polynomials is presented. 
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o In another paper |56| he, with a cohaborator, noted that weak gravita- 
tional, electromagnetic, neutrino and scalar fields, considered as perturbations 
on Kerr background satisfied the Teukolsky Master Equation. Two non-trivial 
equations were obtained after separating the variables, one equation only with 
the polar angle and another using only the radial variable. These were solved 
by transforming each one into the form of a confluent Heun equation. 

Among other papers on this subject one may cite the following papers: 

o In reference [57], R- Manvelyan, H.J.W. Miiller Kirsten, J.Q. Liang and 
Y. Zhang, calculated the absorption rate of a scalar field by a D3 brane in ten 
dimensions in terms of modified Mathieu functions, and obtained the S-matrix 

o In reference |58j, T. Oota and Y. Yasui studied the scalar laplacian on 
a wide class of five dimensional toric Sasaki-Einstein manifolds, ending in two 
Heun's differential equations. 

o In reference [59| . S. Musiri and G. Siopsis found out that the wave equation 
obtained in calculating the asymptotic form of the quasi-normal frequencies for 
large AdS black holes in five dimensions reduced to a Heun equation. 

o A. Al-Badawi and I. Sakalli studied the Dirac equation in the rotating 
Bertotti- Robinson spacetime [60j , ending up with a confluent Heun equation for 
the angular part of the equation for a massless particle. 

o Mirjam Cvetic, and Finn Larsen studied grey body factors and event hori- 
zons for rotating black holes with two rotation parameters and five charges in 
five dimensions. When the Klein Gordon equation for a scalar particle in this 
background is written, one gets a confluent Heun equation. In the asymptotic 
region this equation turns into the hypergeometric form [6lj. When they stud- 
ied the similar problem for the rotating black hole with four U{1) charges, they 
again obtained a confluent Heun equation for the radial component of the Klein 
Gordon equation, which they reduced to the hypergeometric form by making ap- 
proximations . These two papers are partly repeated in [53] . Same equations 
were obtained which were reduced to approximate forms which gave solutions 
in the hypergeometric form. 

Other relevant references I could flnd are listed as references [53] - [TTj . 

I first "encountered" this type of equation when we tried to solve the scalar 
wave equation in the background of the Nutku helicoid instanton 261 ■ this 
case one gets the Mathieu equation which is a special case of the Heun class 

m- 

o The helicoid instanton is a double-centered solution. As remarked above, 
for the simpler instanton solution of Eguchi-Hanson 2i , hypergeometric solu- 
tions are sufficient [25] . Here one must remark that another paper using the 
Eguchi-Hanson metric ends up with the confluent Heun equation |72j . These 
two papers show that sometimes judicious choice of the coordinate system and 
separation ansatz matters. 

o Sucu and Unal obtained closed solutions for the spinor particle written in 
the background of the Nutku helicoid instanton metric p5] . One can show that 
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these solutions can be expanded in terms of Mathieu functions if one attempts 
to use the separation of variables method, as described by L.Chaos-Cador and 
E. Ley-Koo [73]. 

In the subsequent sections I will summarize some work Tolga Birkandan and I 
have done on this topic [Til [75] . 



2 Dirac equation in the background of the Nutku 
helicoid metric |74L 175] 



The Nutku helicoid metric is given as 



ds^ = -^^=[dr^ + ir^+a^)de^+[l + ^sm^e)dy^ 

+ ^ V ^ / 

- ^ sin 29dydz+ ^1 + ^ cos^ 9^ dz^] . (25) 

where Q < r < oo, Q < < 2tt, y and z are along the Killing directions and will 
be taken to be periodic coordinates on a 2-torus |28j . 
If we make the following transformation 



— a sinha;, (26) 



the metric is written as 



a2 

ds'^ = —^\Vih.2x{dx^ ^ dd^ 



+ . [(sinh^ X + sin^ 6>)d;/^ (27) 
smh2x 

- sin 29dydz + (sinh^ x + cos^ 9)dz'^]. 

We write the system in the form L-tp — Aip, where L is the Dirac operator, A is 
the eigenvalue, and try to obtain the solutions for the different components. 
We use the Newman-Penrose formalism [7S1 [77] to write the Dirac equation, 
with four components, in this metric. We use a solution of the type ipi = 
g*(fcyy+fc2 2)v[/j(a;^ 0). We also use ky = kcos(j)^kz = ksin(j). The transformation 
^12 = I ^, fi 2 is used for the upper components to have similar equations 

' y sinhlx ■ 

for all components. Then these equations read: 

(9^+i5e)*3 +m/c[cos(6i- + ia;)]«'4 = A-^/i, (28) 

v2 



(a, - 09)*4 - iak{cos{B - - zx)]*3 - (29) 

v2 
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dStTlh^X 

i-d, + ide)fi +iak[cosi0-(j) + tx)]f2^~A ^^-3, (30) 

v2 



fi mn hyT 

i-d, - ide)f2 - iak[cos{0 - - ix)]h = -A ^^4. (31) 

v2 

We solve our equations in terms of /i^2 and substitute these expressions in 
equations, given above. This substitution gives us second order, but uncoupled 
equations for the lower components: 



dxx + dgg H — 7r[k {—cos [2{e - (j))] - cosh2x) - A^sinh2x]J ^3,4 = 0. (32) 

We can separate this equation into two ordinary differential equations by the 
ansatz 4'3_4 = R{x) S {9 — (f>) . Using this ansatz gives us two ordinary differential 
equations. The equation for S reads 



deeS{e) - (^y fc2cos(2e) - S{e) = 0, 



(33) 



where Q = {9 — (j)). This equation is of the Mathieu type and the solution can 
be written immediately. 

S{9) = CiSe{n, ^,9 - (j)) + C^Soin, ^,9 - 0), (34) 

where Se, So are the even and odd Mathieu functions. The solutions should be 
periodic in the angular variable 8. This fact forces n, the separation constant, 
to take discrete values. It is known that the angular Mathieu functions satisfy 
an orthogonality relation such that functions with different n values are per- 
pendicular to each other. Here, we integrate the angular variable from zero to 
27r. 

The equation for R{x) reads 

dxx - [y (fc2cos/i2x) + A'^sinh2x] + n > R{x) = 0. (35) 

This solution is of the double confluent form which can be reduced to the form 

R{x) = DiSe{n, Aq, i{x + b)) + D2So{n, Ae, i{x + b)). (36) 

by several transformations. Here Ci,C2, Di, D2 are arbitrary constants. A^jb 
are given in the original reference |74j . 

As a result of this analysis we see that the solutions of the Dirac equation, 
written in the background of the Nutku helicoid metric, can be expressed in 
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terms of Mathieu functions, which is a special form of Heun function [731 175] . 
Mathieu function is a related but much more studied function with singularity 
structure sme as the double confluent Heun equation. 

o One can also show |74[ [75] that one can use the similar metric in five 
dimensions and obtain, in general, double confluent Heun functions which can 
be reduced to Mathieu functions. 



3 Scalar field in the background of the extended 
Eguchi-Hanson solution [29] 

To go to five dimensions, we can add a time component to the Eguchi-Hanson 
metric [53] so that we have 

ds' ^ -dt' + -^dr' + r^al + +r'{l- (37) 

where ^ 

o'x = 2 ( ~ cos £^d9 — sin 9 sin S^dcj)) (38) 

ay — i (sin ,^c?6' — sin^cos^d^) (39) 

^^{-d^- cos ed(j)). (40) 

This is a vacuum solution. 

If we take the solution to the wave equation $ as 

^ ^ e'fc*e"'^e'('"+^)«(p(r, 0), (41) 
we can write the reduced wave equation as 

Av?(r, 9) = ( :^drr + ^ c»r + k^r^ + -J 4 + 



%mn cos 9 ~ Mm? + n?) , , , , , 

4dee+Acot9de + '-)(p{r,9) = 0. (42) 

sm 9 

If we take (p{r,9) ~ f{r)g{9), the solution of the radial part is expressed in 
terms of confluent Heun {He) functions. 

/ (r) = i-a' + r^) ^ " Ha (o, m, m, i k'a' , i ™^ - 1 A - J fc^a^ 
+ {a' + r^) " (r^ - a^) ^ " He (o, -m, m, 1 fc^a^ i - i A - 1 fc^a^ 

(43) 
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If the variable transformation r ~ aV cosh x is made, the solution can be ex- 
pressed as 



/ ix) = -^{(siiih i^)r^' He (o, m, m, i Pa^, i - 1 A - J k^a^, 1 cosh (x) + £ 

+ (2 cosh (a;) + 2)^^ ™+3 (2 cosh (.x) - 2) ^ '"+^ (44) 
X He (^0, -TO, TO, i /c^^^ 1 - 1 A - i fc^a^, i cosh (x) + }. 

2 . 2 

We tried to express the equation for the radial part in terms of w = ^^^5— to 
see the singularity structure more clearly. Then the radial differential operator 
reads 

du'^ \u — 1 u J du \u — 1 u J — u)"^ 

This operator has two regular singularities at zero and one, and an irregular 
singularity at infinity, the singularity structure of the confluent Heun equation. 
This is different from the hypergeometric equation, which has regular singular- 
ities at zero, one and infinity. 

The angular solution is in terms of hypergeometric functions. 



1 , /I 

9 (0) = —qW^ - 2 cos (e) i^- cos (0) - -j i^-cos{0) 



[(2 cos (6*) 4-2) "-2™ 

X 2Fi([-n+ivA+l + i, -n-^\/A + l + ^],[l-n-TO],i cos(0) + i) 

(46) 

+ (2 cos (61) + 2)5 + 5 "+5"^ 
X2i^j([m + ivA+l + i, TO- ivA + l + i], [l + n + TO],i cos (61) + i)]}. 



4 Derivation of a new generalized Heun equa- 
tion [29j 

The Dirac equation, with A = 0, written in the background of the Nutku heli- 
coids metric, reads as follows: 

{d^+ide)'i>3 + iak[cos{e - cj) + ix)]-^4 = 0, (47) 



12 



(9:r-«5e)*4 ~iak[cos{e~(j)~ix)]-^3=0, (48) 

{~d^+ide)fi +iak[cos{e-(l) + ix)]f2=0, (49) 

(-a, - idg)f2 - iak[cos{0 - - ix)]fi = 0. (50) 

These equations have simple solutions '251 which can also be expanded in terms 
of products of radial and angular Mathieu functions (73] , [Zl] . Problem arises 
when these solutions are restricted to boundary [^Hl [7S] . 

To impose these boundary conditions we need to write the little Dirac equation, 
the Dirac equation restricted to the boundary, where the variable x takes a fixed 
value xo- We choose to write the equations in the form, 

— {i-j7^^3 +ikacos{e - (l) + ixo)'^4} ^ Xfi, (51) 
a do 

— {"*:^*4 -iakcos{e-(l)-ixQ)^3} = Xf2, (52) 
a da 

— i-'lEf^ -^akco.sie-(j) + ixo)f2} = X'fs, (53) 
a d6 

— {^-JEh +iakcos{0- (t)-ixo)fi} ^ X^4:- (54) 
a dO 

Here A is the eigenvalue of the little Dirac equation. We take A = as the 
simplest case. The transformation 

e = e-^~ixo (55) 
can be used. Then we solve /i in the latter two equations in terms of f2'- 

^/2-tane-/2+ — 



-j^/2-tane — /2+^^[cos(2e)cosh(2a;o)-«sin(2e)sinh(2a;o)+cosh(2a;o)]/2 = 0. 



(56) 

When we make the transformation 

u = e^*®, (57) 

the equation reads, 

{4(u+l)M[u^ + ^]-2zM(u-l)^+.^(ii+l)[ue-2-o + ie2-o+cosh(2xo)]}/2 = 0. 
dw^ du du 2 u 

(58) 

This equation has irregular singularities at m = and oo and a regular singularity 
at u = —1. This is still another generalized Heun equation, different from the 
one given by by reference [46] . 
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